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Abstract 

Non-holonomic deformations of integrable equations of the KdV hierarchy are stud- 
ied by using the expansions over the so-called "squared solutions" (squared eigen- 
functions). Such deformations are equivalent to perturbed models with external 
(self-consistent) sources. In this regard, the KdV6 equation is viewed as a spe- 
cial perturbation of KdV equation. Applying expansions over the symplectic basis 
of squared eigenfunctions, the integrability properties of the KdV hierarchy with 
generic self-consistent sources are analyzed. This allows one to formulate a set 
of conditions on the perturbation terms that preserve the integrability. The per- 
turbation corrections to the scattering data and to the corresponding action-angle 
variables are studied. The analysis shows that although many nontrivial solutions 
of KdV equations with generic self-consistent sources can be obtained by the Inverse 
Scattering Transform (1ST), there are solutions that, in principle, can not be ob- 
tained via 1ST. Examples are considered showing the complete integrability of KdV6 
with perturbations that preserve the eigenvalues time-independent. In another type 
of examples the soliton solutions of the perturbed equations are presented where 
the perturbed eigenvalue depends explicitly on time. Such equations, however in 
general, are not completely integrable. 



1 Introduction 

Nonholonomic deformations of integrable equations attracted the attention of the scientific 
community in the last few years. In [22], based on the Painleve test, applied to a class 
of sixth-order nonlinear wave equations, a list of four equations that pass the test was 
obtained. Among the three known ones, there was a new equation in the list (later known 
as sixth-order KdV equation, or just KdV6): 

+ Vxdx + \ Vx y ( v * + Vxxx ~ 3 ^ = °- W 

Recognising the potential KdV (PKdV) equation in the second factor in ([I]) and the 
Nijenhuis recursion operator for PKdV A/pKdv — ( — \d x + v x d x + \v xx )d~ 1 , after change 
of variables v = v x and w = Vt + v xxx — 3v%, one can convert ([1]) into a "potential" form: 

ut + u xxx - 6uu x -w x = 0, (2) 
-\w xxx + uw x + \wu x = 0, (3) 
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or equivalently 



(^~\ d * +u + \ u $ J ( u t + u xxx - 6uu x ) = 0. (4) 



Here d x 1 is a notation for the left-inverse of d x : d x 1 f(x) = f x f(y) dy. Note also, that d x 1 
is a Hamiltonian operator and — \d x + v x d x + symplectic operator. Since both 

operators —\d x + v x d x + \v xx and A/pKdv are weakly nonlocal, then eqn. (jTJ possesses 
the same recursion operator as the PKdV equation [22| H5] and therefore infinitely many 
integrals of motion. In the same paper [22], a Lax pair and auto-Backlund transformations 
are obtained for KdV6 equation, ([[]). In [13], applying Hirota bilinear method to KdV6 
equation, another class (much simpler) of auto-Backlund transformation is obtained. 
In [33] B. Kupershmidt described (J2J) and (j3J) as a nonholonomic of the KdV equation, 



written in a bi-Hamiltonian form: 

v t = B 1 (^jf) ~ B\w) = B 2 ( 5 -^)- B\w), B\w) = (5) 



where 



B l = d xi B 2 = - l -dl + v x d x + l -v. 



XX 



^ x ■ ' 2 

are the two standard Hamiltonian operators of the KdV hierarchy, and 

i 2 2 3 2 

are the conserved densities for the same hierarchy. The recursion operator for the KdV 
hierarchy is given (after a suitable rescaling) by A = B 2 o (S 1 ) -1 . A more general setup 
suitable for non-Hamiltonian systems of equations is presented in [4"5] . 

Later on, it was shown in [13118] that the KdV6 equation is equivalent to a Rosochatius 
deformation of the KdV equation with self-consistent sources. Soliton equations with 
self-consistent sources have many physical applications [321 HOI ESI HE] • For example, they 
describe the interaction of long and short capillary-gravity waves |39j. Different classes of 
exact solutions of <^ are obtained in [22j [29j |30j [311 [32] . A geometric interpretation of 
the KdV6 equation is given in [17]. 

In the present paper, exploiting the potential form of the KdV6 equation (J2]), we study 
the class of inhomogeneous equations of KdV type 

u t + u xxx - 6uu x = W x [u] (x) , (6) 

with an inhomogeneity/perturbation that presumably belongs to the same class of func- 
tions as the field u(x) (i.e. decreasing fast enough, when \x\ — > oo). 

Generally speaking, the perturbation, as a rule destroys the integrability of the con- 
sidered nonlinear evolution equation (NLEE). The idea of perturbation through non- 
holonomic deformation, however, is to perturb an integrable NLEE with a driving force 
(deforming function), such that under suitable differential constraints on the perturbing 
function(s) the integrability of the entire system is preserved. In the case of local NLEE's, 
having a constraint given through differential relations (not by evolutionary equations) is 
equivalent to a nonholonomic constraint. 

To the best of our knowledge, the most natural and efficient way for studying in- 
homogeneities/perturbations of NLEE integrable by the inverse scattering method is by 
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using the expansions over the so-called "squared solutions" (squared eigenfunctions). The 
squared eigenfunctions of the spectral problem associated to an integrable equation rep- 
resent a complete basis of functions, which helps to describe the Inverse Scattering Trans- 
form (1ST) for the corresponding hierarchy as a Generalized Fourier transform (GFT). 
The Fourier modes for the GFT are the Scattering data. Thus all the fundamental prop- 
erties of an integrable equation such as the integrals of motion, the description of the 
equations of the whole hierarchy and their Hamiltonian structures can be naturally ex- 
pressed making use of the completeness relation for the squared eigenfunctions and the 
properties of the corresponding recursion operator. 

This approach was developed first for the Zakharov-Shabat system [2H [251 [26J [271 
[28| IT3j HI] . In particular, in [13] a special combination of squared solutions named the 
'symplectic basis' was introduced. Its special property consists in mapping the variation 
of the potential of the Zakharov-Shabat system into the variations of the action-angle 
variables. 

Later such basis was introduced also for the Sturm-Liouville's operator [281 E [18] - 
Note that the scattering data of the Sturm-Liouville's operator in a generic situation have 
a pole for k = and their variational derivatives with respect to the potential do not 
vanish when \x\ — > oo. This can lead to to nontrivial contributions to the associated 
symplectic and Poisson structures, coming from the boundary terms [U [9]. For more 
details, we refer to [HI [H] and the references therein. 

The expansions over the squared solutions have been used for the study of the pertur- 
bations of various completely integrable systems, see [261 HTJ [27J . As important particular 
cases these authors analyzed NLEE with singular dispersion relations and have shown 
that whenever the dispersion relation has a pole located on a given eigenvalue A of the 
corresponding Lax operator L, then the corresponding evolution is no longer isospectral. 
Further generalizations are made in [121 E] • An alternative approach to the same type 
of equations has been proposed by Leon in [311 S3 EH [38] ; it is based on the equivalence 
of the Zakharov-Shabat system to a <9-problem. 

The structure of the present paper is as follows: In Section 2 we give some back- 
ground material about the direct scattering problem for KdV hierarchy. In Section 3 we 
describe briefly the generalized Fourier transform for the equations of KdV hierarchy as 
a key point for their integrability. This includes the minimal set of scattering data, their 
time evolution, the expansions over the complete set of squared solutions and over the 
symplectic basis, action-angle variables, etc. In the next Section 4 we formulate general 
conditions on the perturbations (driving forces) for preserving integrability, treating sep- 
arately integrable and non-integrable cases. We provide also a set of nontrivial examples 
for both cases. Then, in Section 5 we present examples for 1-soliton solutions of some of 
the perturbed systems from Section 4. 

2 Preliminaries: Integrability of the Equations of KdV 
Hierarchy 

The spectral problem for the KdV hierarchy is given by the Sturm-Liouville equation 



[42. 



EB] 



(7) 



3 



where u(x) is a real-valued (Schwartz-class) potential on the whole axia and k G C is 
spectral parameter. The continuous spectrum under these conditions corresponds to real 
k. We assume that the discrete spectrum consists of finitely many points k n = in n , 
n = 1, . . . , N where K n is real. 

The direct scattering problem for ((7j) is based on the so-called "Jost solutions" f + (x, k) 
and f + {x, A;), given by their asymptotics: x — » oo for all real k 7^ [42] : 

lim e _ite / ± (a;, k) = 1, fceM\{0}. (8) 

From the reality condition for u(x) it follows that f (x, fc) = f ± (x, —k). 

On the continuous spectrum of the problem ((7j) and the vectors of the two Jost solutions 
are linearly related 0: 

f~(x,k) = a(k)f + (x,-k) + b(k)f + (x,k), Im k = 0. (9) 

In addition, for real k 7^ we have: 

f ± (x,k) = f ± (x,-k). (10) 

The quantities TZ ± {k) = b(±k) / a(k) are known as reflection coefficients (to the right with 
superscript (+) and to the left with superscript (— ) respectively). It is sufficient to know 
TZ ± {k) only on the half line k > 0. Furthermore, TZ ± (k) uniquely determines a(k) [4"2"] . 
At the points k„ of the discrete spectrum, a(k) has simple zeroes i.e.: 

a(k) = (k - iK n )a n + ^(k - m n ) 2 a n H , (11) 

Here, the dot stands for a derivative with respect to k and a n = a(iK n ), a n = a(in n ), etc. 
At the points of the discrete spectrum /" and / + are again linearly dependent: 

f~(x,iK n ) =b n f + (x,iK n ). (12) 

The discrete spectrum is simple, there is only one (real) linearly independent eigenfunc- 
tion, corresponding to each eigenvalue ire n , say f~(x) = f~(x,iK n ). 
The sets 

s ± =\n ± (k) (k>o), n n} Rt= b -^, n = l,...Jv) (13) 

are known as scattering data. Each set - S + or S~ of scattering data uniquely determines 
the other one and also the potential u(x) [32j [181 HH]- 

KdV equation appeared initially as models of the propagation of two- dimensional 
shallow water waves over a flat bottom. More about the physical relevance of the KdV 
equation can be found e.g. in (U [191 EB 171 120]. 



3 Generalized Fourier Transforms 

A key role in the interpretation of the inverse scattering method as a generalized Fourier 
transform plays the so-called 'generating' (recursion) operator: for the KdV hierarchy it 
has the form pQ: 

1 1 r x 

L± = --d 2 + u{x) - - / dxu'{x) ■ . (14) 

1 According to the notations used in [32] f + { x , k) = i>(x, k), f~(x, k) = ip(x, k). 
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The eigenfunctions of the recursion operator are the squared eigenfunctions of the spectral 
problem (J7J): 

F ± (x,k) = k)) 2 , F±(x) = F(x,iK n ), (15) 

From the completeness of the squared eigenfunctions it follows that every function g(x), 
belonging to the same class of functions as the potential u(x) of the Lax operator, can be 
expanded over the two complete sets of "squared solutions" [18] : 



1 



N 



~gHk)F±(x } k)dk =f [9iAU x ) + 92,A 



2ni 

where F^(x) = [■^F ± (x, k)]k=kj an d the Fourier coefficients are 



3,% 



(16) 



9i,j 



ka 2 
1 



1 f°° 
T(Y}(9,F T ), where (g,F) = J g(x)F(x)dx, 



kjhj 



92,j 



kj(l 2 



9,F? 



1 'dj 

kj dj 



Here we assume that in addition / + and /~ are not linearly dependent at x — 0. The 
details of the derivation can be found e.g. in [HI H8] . 

In particular one can expand the potential u(x), the coefficients are given through the 
scattering data [HI EE] : 



u\x) 



±— y kn ± (k)F ± (x, k)dk + 4i 5Z k i R t F t { 



x 



(17) 



The variation fox (a;) under the assumption that the number of the discrete eigenvalues is 
conserved is 



Su(x) = — 



1 f°° N r 

- / 8TZ ± (k)F^(x, k)dk ±2J2\ RpkjF^(x) + 5RfF, 



± 

3,x 



An important subclass of variations are due to the time-evolution of u, i.e. effectively 
we consider a one-parametric family of spectral problems, allowing a dependence of an 
additional parameter t (time). Then 5u(x, t) = u t St + Q((5t) 2 ), etc. The equations of the 
KdV hierarchy can be written as 



u t + d x Q(L±)u(x, t) = 0, 



(19) 



where fl(k 2 ) is a rational function specifying the dispersion law of the equation. The 
substitution of ffl8|) and ffTTj) in ffl9l) gives a system of trivial linear ordinary differential 
equations for the scattering data: 



nf ± 2ikn(k 2 )n ± = o, 

Rf t ± 2ik 3 n{k 2 )Rf = 0, 



\i.t 



0. 



(20) 
(21) 
(22) 



The KdV equation u t — 6uu x + u xxx = can be obtained for Q(k 2 



-4k 2 . 
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Once the scattering data are determined from (1201) - fl22l) one can recover the solution 
from ( jT8l) . Thus the Inverse Scattering Transform can be viewed as a GFT. 

For our purposes, it is more convenient to adopt a special set of "squared solutions" , 
called symplectic basis [T3] ITS] . It has the property that the expansion coefficients of the 
potential u(x) over the symplectic basis are the so-called action-angle variables for the 
corresponding NLEE. 

For the KdV hierarchy, the symplectic basis is given by: 

V{x, k) = T {TZ ± {k)F ± (x, k) - TZ ± (-k)F ± (x, -k)) , (23) 
Q(x, k) = TZ-(k)F~(x, k) + TZ + (k)F + (x, k), (24) 

P n (x) = -R5F£(x), Q n (x) = -±-(r+F-(x) - R-F+(x)). (25) 

Its elements satisfy the following canonical relations: 

[[P(h), Q(k 2 )}} = 8(ki - k 2 ), [[V(h),V(k 2 )]\ = [[Q(h), Q(k 2 )]] = 0, (26) 

[[Pm, Qn}} = Smn, [[P m , P„] = [[Qm, Qn}} = 0, (27) 

(&4 > 0, k 2 > 0) with respect to the skew-symmetric product 

i poo poo 

[[f,g]] = - / (f{x)g x {x) - g(x)f x (x))dx = / f(x)g x (x)dx, (28) 



2 . 

The symplectic basis satisfies the completeness relation [18J: 

{x -y)-e { y- x) _ i rr P{X:k)QM _ Q{Xtk)tPM \M 



N 



J2[ P n(x)Qn(y)-Qn(x)Pn(y)), (29) 



n=l 

where f3(k) = 2ikb(k)b(—k). Notice that the integration over k is from to oo. It follows 
that every function X(x) from the same class as the potential u(x) (i.e. smooth and 
vanishing fast enough when x — > ±00) can be expanded over the symplectic basis: 

1 f°° dk 

X{x) = —J — (V(x, k)Mk) - Q(x, k)p x {k)) 

N 



^ {P n {x)(f>n,X ~ Qn{x)p n ,x) ■ (30) 



n=l 



The expansion coefficients can be recovered from the so-called inversion formulas: 
:(*) = \[Q(y,k),X(y)l, Px (k) = [[V(y,k),X 



X\ / LL — \a i / i w i jj / ' '» v / LL' v " ' >< v " ' JJ ' /o-| \ 

0„,x = [[Q n (y),X(y)}}, p n , x = [[P n (y),X(y)}}. 

In particular, if X(x) = u(x) is a solution of the spectral problem, one can compute [15] : 

[[V(y,k),u(y)]\ = 0, [[Q(y,k),u(y)]\ = -4ik(3(k), (32) 

[[P n (y),u(y)j = 0, [[Qn(y),u(y)j = Mk n . (33) 
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Thus, from fl32l one gets: 



u{ x) 



2 

Tli 



N 







V(x,k)dk -^4ik n P n 



x). 



(34) 



n=l 



The expression for the variation of the potential is 

"°° dk 



5u(x) 



1 
2~ix 



{V x {x,k)5(j>{k)-Q x {x, k)5p(k)) 



N 



^ (Pn,x$4>n ~ Qn,x8p n ) (35) 



n=l 



with expansion coefficients 



p(k) = -2ik\n\a(h)\ = -2ik\n(l-TZ-(k)K-(-k)), k > 0, (36) 

ft-(*o 



0(fc) = 2i/3(k)argb(k) = (3(k) In 



7e+(fc)' 



2 In 6 n = In 



Rf 



(37) 
(38) 



These are known as action (p(/c)) - angle (4>{k)) variables for KdV equation |49jll2]. They 
satisfy the canonical relations 



{ P (k 1 ), ( f ) (k 2 )} = 6(k 1 -k 2 ) 



{ P {k 1 ), P {k 2 )} = {<f)(k 1 ),<f)(k 2 )} = 0, 
{Pm,Pn} = {0m, 0n} = 0. 



(39) 
(40) 



{ki > 0, k 2 > 0) with respect to the Poisson bracket: 

{A,B} = 



SA SB 

5u ' 5u 



One can verify that the scattering data determines the canonical variables and vice- versa. 
The symplectic basis consists of the eigenfunctions of the operator A = \{L + + LJ): 



AV(x, k) = k 2 V(x,k), 
AQ(x,k) = k 2 Q(x,k), 



AP n (x) = KP n (x), 



AQn(x) = k n Q n (x). 



(41) 
(42) 
(43) 
(44) 



Again, for variations due to the time-evolution of u, 5u(x,t) = u t 5t + Q((5t) 2 ), etc. the 
equations of the KdV hierarchy 



u t + d x tt(A)u(x,t) = 0, 



(45) 



with (13~4"1) and (1351) are equivalent to a system of trivial linear ordinary differential equations 
for the canonical variables (which can be considered as scattering data): 



<t> t = 4i{3(k)n(k 2 ), 

<f> nt = iik n Q(kl), 



Pt(k) = 0, 

Pn,t = 0. 



(46) 
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4 Perturbations to the equations of the KdV hierar- 
chy 

4.1 General perturbations 

Let us consider a general perturbation W^fw] to an equation from the KdV hierarchy: 

u t + d x tt(A)u(x,t) = W x [u\. (47) 

The function W^fu] is assumed to belong to the class of admissible potentials for the 
associated spectral problem (EJ) (Schwartz class functions, in our case). 
The expansion of the perturbation over the symplectic basis is: 

1 f°° dk 

W x [u] = — 7T77T (Px(x, k)(j) W {k) - Q x (x, k)p W {k)) (48) 

2ti J fj[k) 

N 

^)0n,iy — Qn,x(%)Pn,W^) • 

n=l 

The substitution of the above expansion (1451) in fj4T|) together with (1541) and (f3"o"j) leads to 



a modification of the time evolution (14*6]) of the scattering data as follows: 

4> t = 4tf3(k)n(k 2 ) + (j) W {k, t- p(k, t), (j)(k, t), Pn (t), 4> n {t)) + <j) Q 5{k), (49) 

pt{k) = pw{k, t; p(k, t), (f)(k, t),p n (t), (j> n (t)), (50) 

(f) n>t = Aik n Vt(k 2 n ) + 4> n ,w(k, t; p(k, t), <p(k, t),p n (t), (f) n (t)), (51) 

p n ,t = pn,w(k, t; p(k, t), (p(k, t), p n (t), 4> n {t)). (52) 



Since W — W[u] and u depend on the scattering data, we observe that the expansion coef- 
ficients of the perturbation {(f> w {k) = \[Q(y, k), W(y)~^ etc.) also depend on the scattering 
data. Thus for generic W the new dynamical system (T4"9"j) - (|5"2"|) for the scattering data 
can be extremely complicated and non-integrable in general. This reflects the obvious 
fact that the perturbed integrable equations are, in general, not integrable. 
Below we consider several important particular cases. 



4.2 Self-consistent sources — integrable case 

It is widely spread that KdV with self-consistent sources is a special case of perturbed 
KdV with 

N 

W x [u] = c n V n , x (x) + c V x (x, 0), (53) 

n=l 

where c n are some constants and the constant Co is usually set to 0. The last term in (|53|) 
affects only the end-point of the continuous spectrum of L. Such choice of the perturbation 
substantially simplifies the evolution of the scattering data: 

fa = 4tP(k)Q(k 2 ) - c S(k), p t (k) = 
0„ >t = Aik n n(k 2 n ) - Cn, Pn,t = 0- 

It is obvious that perturbations (15^1) preserve integrability. Indeed, from eqs. (f5"4"l) one 
finds that such perturbations are isospectral, since k nyt = 0. The only effect of the 
perturbation (1531) consists in changing the time-dependence of the angle variables, related 
to the discrete spectrum of L. 
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Example 1. Letco ^ and all other c n = 0. ThenW x [u] = cqP(x,0). Since AP(x, 0) = 
(see eq. (37\)) we easily conclude, that 

A(v t + v xxx - 3(v 2 ) x - c V(x, 0)) = A(v t + v xxx - 3(v 2 ) x ) = 0, (55) 

i.e such perturbed KdV is equivalent to the KdV6 equation. Thus from the above consid- 
erations we confirm the result in f3^ . that KdV6 is integrable. 

4.3 Self-consistent sources — non-integrable case 

Here we consider a more general type of self-consistent sources given by: 

N 

W x [u\ = ^(c n V n , x (x) + c n Q n>x (x)) + coV x (x, 0), (56) 

71=1 

Again we consider c n , c n and c as constants. Note, that Q x (x, k) is an odd function of k 
and therefore Q x (x, 0) = 0. The evolution of the scattering data under such perturbation 
is: 

(j>t = Ai/3(k)Q(k 2 ) - c S(k), Pt (k) = 

4>n,t = ^ik n Vt{k 2 n ) - Cn, p n>t = C n . 

i.e both the action variables p n and the discrete eigenvalues become time-dependent. Such 
perturbation is not isospectral and obviously violates integrability. In addition the angle 
variables become nonlinear functions of t. 

Let us investigate the integrability of the following equation: 

A*(u t + d x Q{A)u(x, t)) = 0, (58) 

where the star is a notation for a Hermitian conjugation. KdV6 in (j4j) is a particular case 
of this equation with Q(A) = — 4A. In order to simplify our further analysis, instead of 
the equation (|58|) we study the following one: 

(A* -\ 1 )(u t + d x Q(A)u(x,t)) = 0, (59) 

where Ai is a constant. The corresponding analogue for KdV6 is 

v 6x + v txxx - 2v t v xx - iv x v xt - I0v x v 4x - 20v xx v xxx + 30vlv xx 

+ 4Ai(urt + v xxxx - Qv x v xx ) = 0. 

Due to the identity dA = A*d we can represent fl59|) in the form 

8(A - X^dd-'ut) + n(A)u(x,t)) = 0. (60) 

Since the operator d does not have a kernel when u is Schwartz class, (!60|) is equivalent 
to 

ut + d x £l(A)u(x, t) 

where 

Cl 
Cl 



(ciPi(a;,*) + CiQi(a;,*)) a . for \ x = k\ < 0, 

{c 1 V(x,k 1 ,t) + c 1 Q(x,k 1 ,t)) x for Ai = ^> , 



(61) 



Ci(t, fci; p{h,t), <f>(kt, t),p n (t), 4> n (t)), 

cx(t, ki, p(h, t), <f>(k u t),p n (t), (j) n (t)) 
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are ^-independent functions, but the important observation is that the time-dependence 
could be implicit through the scattering data of the potential u(x, t). Equation ( |6T|) is a 
perturbed equation from the KdV hierarchy. The perturbation in the right-hand side of 
(IBTj) is in the eigenspace of the recursion operator corresponding to the eigenvalue Aj., i.e. 
it is given by 'squared' eigenfunctions of the spectral problem ([7]) at Ai. 

If Ai < is a discrete eigenvalue, the corresponding dynamical system (j4"9"|) - (1521) for 
the scattering data is 

<j> t = 4iP(k){l(k 2 ), (62) 

Pt(k) = 0, (63) 

<f> n ,t = 4ik n Q(kl) - ci{t,fa\p{fa,t),<j)(ki,t),p n (t),<f> n (t))5n,i, (64) 

Pn,t = c 1 (t,k 1 ;p(k 1 ,t),(f)(k 1 ,t),p n ( y t),(j) n (t))5 ntl . (65) 



If Ai > is a continuous spectrum eigenvalue, the dynamical system (H9l) - (1521) has the 
form 



h t = 4i/3(k)n(k 2 ) 



+ 2nP(k)c 1 (t, fa; p(fa,t), <j)(fa, t),p n (t), Mt))S(k - fa), (66) 

Pt (k) = -27T/3(A;)ci(t, fci; p(fa,t), <j>{k u t),p n (t), <p n (t))6(k - fa), (67) 

<i> n ,t = tik n Sl{kl), (68) 

Pn,t = 0. (69) 

It is clear, that dynamical systems like ( 1621) - ( !65l) or (166]) - ( l69l can not be integrable for 
a general functional dependence of c\ and ci on the scattering data. Thus the generic per- 
turbed equations from KdV Hierarchy (158 1) , are not completely integrable. In other words, 
there are solutions, which can not be obtained via the Inverse Scattering Method, since 
the aforementioned dynamical systems for the scattering data are not always integrable. 

Example 2. Let us assume that Ai = — k\ < is a discrete eigenvalue of L, and that c\ 
and c\ in eq. ( 05)) are constants. Then the corresponding perturbed KdV will be equivalent 
to the following equations for the scattering data of L: 



t = 4iP(k)n(k z ), Pt(k) = o, 



4dk n Q(k n ) - ciS n>1 , p n>t = ci<5 n)1 . 



(70) 



From the last of these equations we find 

w 2(t) = c 1 t + «?(0), (71) 

i.e the eigenvalue Ki(t) explicitly depends on t. Clearly (71\ ) makes sense for those t for 
which C\t + «f(0) > 0. Then (i) if c\t + «f(0) < the eigenvalue overlaps with the 
continuous spectrum of L and the solution develops singularities; (ii) if for some t K\ 
crosses another discrete eigenvalue the spectrum cease to be simple and the solution is not 
in the assumed Schwartz class, i.e. again may develop singularities. 

Example 3. Let us assume that \\ = —nf < is the discrete eigenvalue of L, and that 
Ci is constant and C\ is a function of time, for example: 

Ci(t) = e K 1 (0)o;cos(a;t) (72) 
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where eo < 1 and uq are constants. Then the eigenvalue k\ depends explicitly on t 
according to 

K \(t) = «?(0)(1 + e sm(ut)). (73) 

// e is chosen small enough the corresponding eigenvalue will not overlap with the con- 
tinuous spectrum, nor with other discrete eigenvalues. Otherwise - as above and again we 
may have singularities. 



5 Examples of exact solutions 

As we have seen, KdV6 in general is not completely integrable system. However, there 
are many solutions of fl58|) or (|60|) that can be written explicitly (apart from the KdV 
solutions, which are an obvious subclass of solutions). 

One such example is when c\ and c\ depend only on t. Then the system ( 162|) - ( |65i) is 
integrable, with new canonical variables 



<fin = 0n(O) + Mkr$l(k n )t - 5 n , :1 / Ci(t)(2t, 

Pn = Pn{0) + S n>1 / ci{r)dr. 



Since c\(t) and c\ do not depend on the scattering data and therefore on u(x,t), the 
Poisson brackets between <p n and p n are canonical. 

This is most often the situation that many authors assume about KdV6 (i.e. c\ = ci(t) 
and Ci = Ci(t) simply given functions of t) when discuss its integrability. 

There are, however other exactly solvable cases. Indeed, if 

TV 

Cl = ^(«lfe0fc + PlkPk), 
k=l 
N 

ci = ^(a 2 fe0fc + fokPk) 

k=l 

are linear combinations of the scattering data with some constants (5 Ptq the system 
(1621) - ( 1651) becomes a system of linear ordinary ODEs and therefore is integrable. 

To illustrate the effect of non-holonomic deformations on the soliton solutions, let 
us consider the 1-soliton solutions corresponding to the specific choice of perturbation 
functions C\ and c\ considered in examples 2 and 3 in the previous Section. The one- 
soliton solution for the KdV equation (corresponding to a discrete eigenvalue Ai = — k\) 
is given by [4"2"] : 

2 k 2 

u ls {x, t) = 7 1 - x , (74) 



ch 2 m (x - AKft - g 



where 0i is the angle variable that corresponds to the prescribed eigenvalue Ai = — K.f. 
Here, also, the quantity v = <Ak\ can be interpreted as a velocity of the soliton, (pi - as a 
position of soliton's mass-center at t — (it is often called a soliton phase). 
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For the special choice of constants in Example 2 one can write down the corresponding 
perturbed 1-soliton solutionin the following form: 



u ls (x,t) 



2(cit + «i(0) : 



Ch 2 \ y/cit+K^Oy 



X 



4k 2 - ^ )t - £i 

1 Kl J Kl 



(75) 



Here the correction in the velocity of the soliton v = v + ^ is due to shift in the angle 
variable 4>i- 

For the case of Example 3 the corresponding perturbed 1-soliton solution reads: 



u ls (x,t) 



2«i(0) 2 (l + e smut) 



ch 2 \ ki(0)V1 + e sinu;t x + 4k 2 (0)(1 + e sinwt)* 



(76) 



Here the soliton velocity depends on time, which is typical for the so-called boomerons 
[3]. Note however, that the boomeron type behavior is characteristic for mult i- component 
KdV equations, while here we get it as an effect of the perturbation. Moreover, the 
time dependence of the velocity is controlled by the time- depending coefficient in the 
perturbation c\{t). 

There are of course many other nontrivial integrable examples that one can construct 
by various choices of the functions c\ and c\ - when ( 162]) - ( 165]) is an integrable system of 
ODEs. In general, however this system is not integrable and therefore there exist solutions 
that can not be constructed by the means of 1ST. 



6 Discussions and Conclusions 

Here we have used the expansion over the eigenfunctions of the recursion operator for the 
KdV hierarchy for studying nonholonomic deformations of the corresponding NLEE from 
the hierarchy. We have shown, that in the case of self-consistent sources, the corresponding 
perturbed NLEE is integrable, but not completely integrable. 

The perturbation results for the Zakharov-Shabat (ZS) type spectral problems have 
been obtained firstly in [23] and for KdV in [23J . As it has been explained, the perturbation 
theory is based on the completeness relations for the squared eigenfunctions. For the 
Sturm-Liouville spectral problem such relations apparently have been studied as early as 
in 1946 [2j and then by other authors, e.g. [2E1 HE]- The completeness relation for the 
eigenfunctions of the ZS spectral problem is derived in [25J and generalisations are studied 
further in [10] [12] [T3] [11] [15] [33], see also [13] (and the references therein). 

The approach presented in this article can be applied also to the study of inhomoge- 
neous versions of NLEE, related to other linear spectral problems, e.g. the Camassa-Holm 
equation [5] E], various difference and matrix generalizations of KdV-like and Zakharov - 
Shabat spectral problems, various non-Hamiltonian systems [4"5] . etc. 
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